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Abstract. The unitary isobar model MAID has been extended and used for a partial wave analysis of 
pion photo- and electroproduction in the resonance region W < 2 GeV. Older data from the world data 
base and more recent experimental results from Mainz, Bates, Bonn and JLab for Q 2 up to 4.0 (GeV/c) 2 
have been analyzed and the Q 2 dependence of the helicity amplitudes have been extracted for a series of 
four star resonances. We compare single-Q 2 analyses with a superglobal fit in a new parametrization of 
Maid2003 together with predictions of the hypercentral constituent quark model. As a result we find that 
the helicity amplitudes and transition form factors of constituent quark models should be compared with 
the analysis of bare resonances, where the pion cloud contributions have been subtracted. 
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1 Introduction 

Our knowledge of nucleon resonances is mostly given by 
elastic pion nucleon scattering [T] . All resonances that are 
given in the Particle Data Tables |2j have been identified 
in partial wave analyses of nN scattering both with Breit- 
Wigner analyses and with speed-plot techniques. From 
these analyses we know very well the masses, widths and 
the branching ratios into the irN and irirN channels. These 
are reliable parameters for all resonances in the 3 and 4 
star categories. There remain some doubts for the two 
prominent resonances, the Roper Pn(1440) which appears 
unusually broad and the <Sii(1535) that cannot uniquely 
be determined in the speed-plot due to its position close to 
the r/N threshold. Both resonances, however, have recently 
been found on the lattice in a very precise calculation with 
a pion mass as low as 180 MeV and converge very close 
to the empirical masses [3]. This has been achieved in a 
quenched calculation giving rise to the conclusion that 
both resonance states are simply qqq states. 

Starting from these firm grounds, using pion photo- 
and electroproduction we can determine the electromag- 
netic -fNN* couplings. They can be given in terms of elec- 
tric, magnetic and charge transition form factors G* E (Q 2 ), 
G* M (Q 2 ) and G* C (Q 2 ) or by linear combinations as he- 
licity amplitudes A 1 / 2 (Q 2 ) 1 A 3 / 2 {Q 2 ) and S 1 / 2 (Q 2 )- So 
far, we have some reasonable knowledge of the A^i% and 
A3/2 amplitudes at Q 2 — 0, which are tabulated in the 
Particle Data Tables. For finite Q 2 the information found 
in the literature is very scarce and practically does not 
exist at all for the longitudinal amplitudes Si/2- But 
even for the transverse amplitudes only few results are 



firm, these are the G* M form factor of the Z\(1232) up to 
Q 2 « 10 GeV 2 , the A 1/2 {Q 2 ) of the 511(1535) resonance 
up to Q 2 w 5 GeV 2 and the asymmetry A(Q 2 ) = (A 2 
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)/{A\ /2 +Al /2 ) for the D 13 (1520) and Fi 5 (1680) reso- 



nance excitation up to Q 2 « 3 GeV 2 which change rapidly 
between -1 and +1 at small Q 2 w 0.5 GeV 2 gj. Fre- 
quently also data points for other resonance amplitudes, 
e.g. for the Roper are shown together with quark model 
calculations but they are not very reliable. Their statis- 
tical errors are often quite large but in most cases the 
model dependence is as large as the absolute value of the 
data points. In this context it is worth noting that also the 
word 'data point' is somewhat misleading because these 
photon couplings and amplitudes cannot be measured di- 
rectly but can only be derived in a partial wave analysis. 
Only in the case of the Z\(1232) resonance this can and 
has been done directly in the experiment by Beck et al. 
at Mainz For the Delta it becomes possible due to 
two important theoretical facts, the Watson theorem and 
the well confirmed validity of the s+p — wave truncation. 
Within this assumption a complete experiment was done 
with polarized photons and with the measurement of both 
7T° and 7r + in the final state, allowing also for an isospin 
separation. For other resonances neither the theoretical 
constraints are still valid nor are we any close to a com- 
plete experiment. The old data base was rather limited 
with large error bars and no data with either target or re- 
coil polarization was available. Even now we do not have 
many data points with double polarization, however, the 
situation for unpolarized e + p — > e' + p + ir° has consid- 
erably improved, mainly by the new JLab experiments in 
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all three halls A,B and C. These data cover a large energy 
range from the Delta up to the third resonance region with 
a wide angular range in 9 V . Due to the 2tt coverage in the 
4> angle a separation of the unpolarized cross section 
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in three parts becomes possible and is very helpful for the 
partial wave analysis. Even without a Rosenbluth separa- 
tion of doT and da^ we have an enhanced sensitivity of 
the longitudinal amplitudes due to the g?ctlt interference 
term. Such data are the basis of our new partial wave anal- 
ysis with an improved version of the Mainz unitary isobar 
model MAID. 



2 Photo- and electroproduction 

For our analysis of pion electroproduction we will use the 
dynamical model DMT [B] and the unitary isobar model 
MAID [Jj . In the dynamical approach to pion photo- and 
electroproduction [8], the t-matrix is expressed as 



t in (E) = u 77r + v J7r g {E) Un(E) . 



(2) 



where v J7r is the transition potential operator for j*N —>■ 
ttN, and t^N and go denote the ttN t-matrix and free 
propagator, respectively, with E = W the total energy in 
the CM frame. A multipolc decomposition of Eq. © gives 
the physical amplitude in channel a [B], 



t<tf(q E ,k;E + ie) =cxp(iS^) cos<J< Q > x [v$(q E) k) 
l q' 2 R^(q E ,q J ;E)v^(q',k). 
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where and R^n are the ttN scattering phase shift and 
reaction matrix in channel a, respectively; q E is the pion 
on-shell momentum and k = |k| is the photon momentum. 
The multipolc amplitude in Eq. © manifestly satisfies the 
Watson theorem and shows that the jtt multipoles depend 
on the half-off-shell behavior of the ttN interaction. 

In a resonant channel the transition potential v 77r con- 
sists of two terms 



v lv (E) 



(4) 



where v?„ is the background transition potential and 
v R 7V (E) corresponds to the contribution of the bare res- 
onance excitation. The resulting t-matrix can be decom- 
posed into two terms jS] 



t^(E) = t*(E)+t*(E), 



where 
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■v^g (E)U N (E), 
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Here i£^. includes the contributions from the nonresonant 
background and renormalization of the vertex j*NR. The 
advantage of such a decomposition is that all the processes 
which start with the excitation of a bare resonance are 
summed up in t R n . Note that the multipole decomposition 
of both t? v and would take the same form as Eq. J3J . 

As in MAID 7 , the background potential v?f(W, Q 2 ) 
was described by Born terms obtained with an energy de- 
pendent mixing of pseudovector-pseudoscalar nNN cou- 
pling and t-channel vector meson exchanges. The mixing 
parameters and coupling constants were determined from 
an analysis of nonresonant multipoles in the appropriate 
energy regions. In the new version of MAID, the S, P, D 
and F waves of the background contributions are unita- 
rized in accordance with the K-matrix approximation, 



t£^ Q (MAID) = cxp (id^) cos 5 {a) v®f(W,Q 2 



(8) 



From Eqs. © and (jHJl, one finds that the difference 
between the background terms of MAID and of the dy- 
namical model is that off-shell rescattering contributions 
(principal value integral) are not included in MAID. To 
take account of the inelastic effects at the higher energies, 
we replace exp (id^) cos 5^ = ^(exp (2iS^) + l) in Eqs. 
© and © by \(r] a exp (2i<5 (a) ) + 1), where r\ a is the in- 
elasticity. In our actual calculations, both the irN phase 
shifts S( a ) and inelasticity parameters r\ a are taken from 
the analysis of the GWU group JU| • 

Following Ref. 0, we assume a Breit-Wigner form for 
the resonance contribution A R (W, Q 2 ) to the total multi- 
pole amplitude, 



2 ,Mw)r R M R u R (w) 
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where f nR is the usual Breit-Wigner factor describing the 
decay of a resonance R with total width r R (W) and phys- 
ical mass M R . The expressions for f jR , f nR and r R are 
given in Ref. 7 . The phase <f>(W) in Eq. © is introduced 
to adjust the phase of the total multipole to equal the 
corresponding nN phase shift £W. While in the original 
version of MAID 7 only the 7 most important nucleon 
resonances were included with mostly only transverse e.m. 
couplings, in our new version all four star resonances below 
W = 2 GeV are included. These are P 33 (1232), P U (144Q), 
£>i 3 (1520), 5 n (1535), P 33 (1232), 5 3 i(1620), S u (1650), 
D 15 (1675), F 15 (1680), £ 33 (1700), Pi 3 (1720), F 35 (1905), 
P 3 i(1910) and P 37 (1950). 

The resonance couplings A^(Q 2 ) are independent of 
the total energy and depend only on Q 2 . They can be 
taken as constants in a single-Q 2 analysis, e.g. in photo- 
production, where Q 2 = but also at any fixed Q 2 , where 
enough data with W and 9 variation is available. Alterna- 
tively they can also be parametrized as functions of Q 2 in 
an ansatz like 



A a (Q 2 ) = A*(0) 



1 + cfQ 2 



(1 



(10) 



with n > 2. Also other parameterizations with an asymp- 
totic fall-off for large Q 2 , e.g. of Gaussian form work 
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Table 1. Recent experimental data of ir electroproduction on 
the proton. The Mainz experiment was done with beam and 
recoil polarization, all others are unpolarized measurements. 
From JLab Hall B data sets at fixed Q 2 of 0.4, 0.525, 0.65, 
0.75 0.90, 1.15 and 1.45 GeV 2 have been used. 



laboratory 


Q\GeV 2 ) 


w cm 


(MeV) 


0T\deg) 


Mainz 20 
Bates [2T] 


0.121 


1232 




180° 


0.126 


1152 


- 1322 


0-38° 


Bonn [22] 


0.630 


1153 


- 1312 


5 - 175° 


JLab, Hall A|2l] 


1.0 


1110 


- 1950 


146 - 167° 


JLab, Hall B|2l) 


0.4-1.8 


1100 


- 1680 


26 - 154° 


JLab, Hall Cl231 


2.8,4.0 


1115 


- 1385 


25 - 155° 



4+ = 5V /p)(% - M t +) , 
S$ = J ^ l+ (14) 
and for j = (£ + 1) - 1/2 

4/t 1)_ = \W + 2)M(^+i)- " , 

Mt+1)- _ £+ 1 e 

b i/2 l=-Z>(t + i)-. (15) 



equally well and can effectively be rewritten in either form. 
With such an ansatz it is possible to determine the param- 
eters A a (0) from a fit to the world database of photopro- 
duction, while the parameters cf and can be obtained 
from a combined fitting of all electroproduction data at 
different Q 2 . The latter procedure we call the 'superglobal 
fit'. In MAID the photon couplings A a are direct input 
parameters. 

Eq.(9) can also be used for a general definition. At the 
resonance position, W = Mr we obtain 

A* (M R ,Q 2 ) = iA a (Q 2 )f lR (M R )U R (MB.)c 7vN e^ M « ) 
= A r a es (M R ,Q 2 )e^ MR 1 (11) 



with f lR {M R ) = 1 and 



1 k w m N r^N 
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The factor c^n is \/3/2 and — 1/\/3 for the isospin 3/2 
and isospin 1/2 multipoles, respectively. This leads to the 
definition 



A a {Q 2 ) 



1 



■lmA r a es (M R ,Q 2 



(13) 



It is important to note that by this definition the phase 
factor e 1 ^ in Eqs. (9 and 11) is not considered as part of 
the resonant amplitude but rather as an artifact of the 
unitarization procedure. In the case of the zi(1232) reso- 
nance this phase vanishes at the resonance position due 
to Watson's theorem, however, for all other resonances it 
is finite and in some extreme cases it can reach values of 
about 60°. A a is a short-hand notation for the electric, 
magnetic and longitudinal multipole photon couplings of 
a given partial wave a. As an example, for the P33 partial 
wave the specific couplings are denoted by E 1+ , M 1+ and 
S\ + . By linear combinations they are connected with the 
more commonly used helicity photon couplings ^1/2, ^3/2 
and Si/2- 

For resonances with total spin j = £ + 1/2 we get 



A l/2 



1 



[{£ + 2)E i+ + £M i+ ] 



With some care, as will be discussed below, the photon 
couplings can be directly compared to matrix elements of 
the electromagnetic current calculated in quark models 
between the nucleon and the excited resonance states, 



/2 
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where J + = 



-^(Jx +iJy) ■ However, these couplings are 
only defined up to a phase £. Since the sign of the pionic 
decay of the resonance has been ignored in the empirical 
definition of these amplitudes, Eq. (13), it must be taken 
into account in a model calculation in order to make com- 
parison with the empirical data. Therefore the phases £ 
have to be individually calculated in each model. In the 
calculations found in the literature this has often been 
ignored and causes some confusion in comparing these 
numbers, especially in critical cases as for the Roper reso- 
nance Pn(1440), where the correct sign cannot simply be 
guessed. 



3 The hypercentral constituent quark model 

As an application for evaluating the photon couplings 
we have used the hypercentral Constituent Quark Model 
(hCQM) JT]. It consists of a hypercentral quark interac- 
tion containing a linear plus coulomb-like term, as sug- 
gested by lattice QCD calculations 



V(x) = - 



with x = y p 



(17) 



where x is the hyperradius defined in terms of the stan- 
dard Jacobi coordinates p and A. We can think of this 
potential both as a two-body potential in the hypercentral 
approximation or as a true three-body potential. A hyper- 
fine term of the standard form J3| is added and treated 
as a perturbation. The parameters a, r and the strength 
of the hyperfine interaction are fitted to the spectrum 
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(a = 1.61 fm~ 2 , t = 4.59 and the strength of the hy- 
perfine interaction is determined by the A - Nucleon mass 
difference). Recently, isospin dependent terms have been 
introduced ^3 m the hCQM hamiltonian. The complete 
interaction used is given by 



Hi, 



V(x) +H S +H 1 + H S j . 



(18) 



Having fixed the values of all parameters, the resulting 
wave functions have been used for the calculation of the 
photocouplings ^j, the transition form factors for the 
negative parity resonances the elastic form factors 
|17II18| and now also for the longitudinal and tranverse 
transition form factors for all the 3- and 4-star and the 
missing resonances. 
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Fig. 1. The Q 2 dependence of the N — > A helicity amplitudes. 
The solid and dashed curves are the results of the superglobal 
fit with MAID and the predictions of the hyperspherical con- 
stituent quark model. The dotted lines show the pion cloud 
contributions calculated with DMT. The data points at finite 
Q 2 are the results of our single-Q 2 fits, see Table 1 for refer- 
ences. At Q 2 = for A1/2 and A 3 / 2 the photon couplings from 
PDG are shown gj. 



4 Data analysis 

The unitary isobar model MAID was used to analyze 
the world data base of pion photoproduction and recent 
differential cross section data on p(e, e'p)ir° from Mainz, 
Bates, Bonn and JLab. These data cover a Q 2 range from 
0.1 • ■ - 4.0 (GeV/c) 2 and an energy range 1.1 < W < 2.0 
GeV, see Table 1. In a first attempt we have fitted each 
data set at a constant Q 2 value separately. This is similar 
to a partial wave analysis of pion photoproduction and 
only requires additional longitudinal couplings for all the 
resonances. The Q 2 evolution of the background, Born 
terms and vector meson exchange, is described with a 
standard dipole form factor. In a second attempt we have 
introduced a Q 2 evolution of the transition form factors 
of the nucleon to N* and A resonances and have parame- 
terized each of the transverse (A 1 ^ 2 and A 3 / 2 ) an d longi- 
tudinal (Si/ 2 ) helicity amplitudes. In a combined fit with 
all electroproduction data from the world data base of 
GWU/SAID [El and the data of our single-Q 2 fit we ob- 
tained a Q 2 dependent solution (superglobal fit). In Fig. 1 
we show our results for the Z\(1232) excitation. Our super- 
global fit agrees very well with our single-Q 2 fits, except 
for the 2 lowest points of S1/2 from our analysis of the 
Hall B data. Whether this is an indication for a differ- 
ent Q 2 dependence has still to be investigated. Generally, 
all our single-Q 2 points are shown with statistical errors 
from x 2 minimization only. A much bigger error has to be 
considered for model dependence. 

We also compare our empirical analyses with the pre- 
dictions of the hypercentral constituent quark. It turns 
out that the transverse amplitudes of the quark model are 
about half of the magnitudes and for the longitudinal am- 
plitude Si/2(Q 2 ) the quark model gives essentially zero. 
This is due to the fact that in the empirical analysis the 
resonance contribution is fully dressed by the pion cloud 
which is not the case in a constituent quark model. As de- 
picted in Fig. 2 a fully dressed resonance contribution is 
renormalized on each vertex and in the propagator. The 





Fig. 2. Resonances with dressed and bare electromagnetic ver- 
tices. 



baryonic states of the hCQM including hyperfine interac- 
tion can be considered as resonances dressed by hadronic 
interaction giving rise to the empirical masses. However, 
the electromagnetic vertex correction (third part of Fig. 2) 
is not included and has to be calculated separately. We 
have already started to do this and in a first attempt we 
have used the dynamical model DMT and extracted the 
pion loop contributions for the s- and p- waves. This esti- 
mate of the pion cloud vertex correction is shown as dotted 
lines in the figures. For the longitudinal Delta excitation 
the entire amplitude is practically given by the pion cloud 
contribution and only a negligible part arises from a bare 
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Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 

Fig. 3. The Q 2 dependence of the transverse and longitudinal helicity amplitudes for the Pn (1440) and the 5n(1535) resonance 
excitation. The notation of the curves and the data is the same as in Fig. 1. 



Delta. The same is true for the electric amplitude which is 
given by the combination Ai/ 2 — A 3 / 2 / \/3. This calculation 
also explains why none of the constituent quark models 
was ever able to give the empirical strength of the Ml ex- 
citation (or the transition moment una) of the Z\(1232). 
While in a simple SU(2) calculation the transition mo- 
ment is lower by about 30% in more refined calculations 
it can be as low as only half of the empirical value. This 
is also the case here in our hCQM, even if it has more 
realistic wave functions. 

In Fig. 3 we show our results for the helicity amplitudes 
of the Roper resonance Pn(1440) and the Sn (1535). For 
these resonance amplitudes the pion cloud contributions 
are most important near the photon point and become al- 
ready negligible around Q 2 = 0.5 GeV 2 . The comparison 
between the hCQM and the empirical amplitudes is rea- 
sonably good, except for the Ai/ 2 amplitude of the Roper. 
This finding has to be further investigated both in the 
framework of the quark model and also in the empirical 
analysis. Certainly, for the Roper resonance the existing 
data is not very sensitive to this partial wave. Further ex- 
periments with double polarization could be very helpful 
to solve this problem. 

Finally, in Figs. 4 and 5 we show our results for the 
Z?i3(1520) and the F 15 (1680) resonances. Here the largest 
discrepancies between our quark model calculations and 
the empirical analysis appear in the helicity 3/2 ampli- 
tudes at small Q 2 . So far, the dynamical model calcula- 
tions have only been done for s- and p- waves, therefore 
we cannot give a pion cloud calculations for these par- 
tial waves. However, our findings encourages very strongly 
such extensions of the dynamical model. Furthermore we 
also have some empirical results for the partial waves that 
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Fig. 4. The Q 2 dependence of the transverse helicity ampli- 
tudes for the Di3(1520) resonance excitation. The notation of 
the curves and the data is the same as in Fig. 1. 



are not shown here, but most of them come out of the 
fit with rather large errors bars in the single-Q 2 analysis. 
This gives us less confidence also for our superglobal fit. 
The reason for it is mainly that we have fewer data points 
to analyze at higher energies. 
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5 Conclusions 

Using the world data base of pion photo- and electropro- 
duction and recent data from Mainz, Bonn, Bates and 
JLab we have made a first attempt to extract all longitu- 
dinal and transverse helicity amplitudes of nucleon res- 
onance excitation for four star resonances below W — 
2 GeV. For this purpose we have extended our unitary 
isobar model MAID and have parametrized the Q 2 depen- 
dence of the transition amplitudes. Comparisons between 
single-Q 2 fits and a Q 2 dependent superglobal fit give us 
confidence in the determination of the Delta amplitudes. 
We can also reasonably well determine the amplitudes 
of the Pu (1440), 5 U (1535), Z>i 3 (1520) and the P 15 (1680), 
even though the model uncertainty of these amplitudes 
can be as large as 50% for the longitudinal amplitudes 
of the Z?i3 and F15. For other resonances the situation 
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Fig. 5. The Q dependence of the transverse helicity ampli- 
tudes for the Fi5(1680) resonance excitation. The notation of 
the curves and the data is the same as in Fig. 1. 



is even worse. However, this only reflects the fact that 
precise data in a large kinematical range are absolutely 
necessary. In some cases double polarization experiments 
are very helpful as has already been shown in pion pho- 
toproduction. Furthermore, without charged pion electro- 
production, some ambiguities between partial waves that 
differ only in isospin as 6*11 and S31 cannot be resolved 
without additional assumptions. Finally, all results dis- 
cussed here are only for the proton target. We have also 
started an analysis for the neutron, where much less data 
are available from the world data base and no new data 
has been analyzed in recent years. Since we can very well 



rely on isospin symmetry, only the electromagnetic cou- 
plings of the neutron resonances with isospin 1/2 have to 
be determined. We have found such a solution for the neu- 
tron and will implement it in the next version of MAID 
(MAID2003). It will be a challenge for the experiment to 
investigate also the neutron resonances in the near future. 
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